Thermal Phase Transitions and Gapless Quark Spectra 
in Quark Matter at High Density 
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Thermal color superconducting phase transitions in three-flavor quark matter at high baryon 
density are investigated in the Ginzburg-Landau (GL) approach. We constructed the GL potential 
near the boundary with a normal phase by taking into account nonzero quark masses, electric 
charge neutrality, and color charge neutrality. We found that the density of states averaged over 
paired quarks plays a crucial role in determining the phases near the boundary. By performing 
a weak coupling calculation of the parameters characterizing the GL potential terms of second 
order in the pairing gap, we show that three successive second-order phase transitions take place 
as the temperature increases: a modified color-flavor locked phase {ud, ds, and us pairings) — » a 
"dSC" phase (ud and ds pairings) — > an isoscalar pairing phase (ud pairing) — * a normal phase 
(no pairing). The Meissner masses of the gluons and the number of gapless quark modes are also 
studied analytically in each of these phases. 

PACS numbers: 12.38.-t,12.38.Mh,26.60.+c 



I. INTRODUCTION 

Over the past few years, the properties of color su- 
perconducting quark matter and the phase structure at 
high baryon density, which were originally studied in ear- 
lier publications Q|, have been examined intensively (see, 
e.g., Ref. 0] for reviews). A color superconductor pre- 
dicted to occur in weak coupling is a relativistic system 
in which the long-range color magnetic interaction, which 
is screened only dynamically by Landau-damping of ex- 
changed gluons, is responsible for the formation of the 
superconducting gap. The gap in weak coupling has a 
non-BCS form, A cx fi exp(— Sir 2 /V2g) with g the strong 
coupling constant and u the quark chemical potential 
Q. Furthermore, the gap has a matrix structure with 
(N c x Nf) 2 components due to various combinations of 
iV c color and Nf flavor degrees of freedom. 

There is a good deal of evidence, obtained from various 
weak coupling analyses 0, IM B ; that if all the quark 
masses are zero (m u — nid = m s = 0) and hence the 
quark chemical potentials are equal (fj, u = Ud — Ms),the 
quark matter is in a color-flavor locked (CFL) phase y| at 
low temperature. A transition from the CFL phase to the 
normal phase in mean-field theory is of second order. In 
weak coupling, the transition temperature T c is related to 
the zero temperature gap At=o as T c — 2 1 / 3 x 0.57At=o 
0, which is a BCS relation except for a factor 2 1 / 3 . 

The massless situation thus described may be approx- 
imately realized at asymptotically high densities where 
the quark masses are negligible compared to the quark 
chemical potentials. However, the effect of nonzero quark 
masses becomes important when the chemical potentials 
decrease. In the presence of a quark mass difference char- 
acterized by 2m s /(m u + md) ~ 25, (3 equilibrium, electric 
neutrality, and color neutrality combine to produce non- 
trivial chemical potential differences between flavors and 



and between colors, 



leading to new phases such as the 
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gapless CFL (gCFL) phase [ly in which two out of nine 
quark quasiparticles are gapless. 

In our recent Letter [ijj , we have studied what kind of 
phase structure appears in /3 equilibrated neutral quark 
matter near the super-normal phase boundary when 
there are quark mass differences. An advantage of study- 
ing the region near the phase boundary is that in clas- 
sifying the possible phase structures we can make use of 
the model-independent Ginzburg-Landau (GL) analysis 
0, IT^ in which the thermodynamic potential difference 
between the superfluid and normal phases is expanded 
in terms of the order parameter (the pairing gap) . Fur- 
thermore, by calculating the parameters controlling the 
thermodynamic potential terms in weak coupling, one 
can prove which phase is realized below the super-normal 
boundary in the high density regime where the weak cou- 
pling analysis is valid. 

A crucial observation found in Ref. is that the 
phase structure near the critical temperature is essen- 
tially dictated by the average density of states of differ- 
ent quarks. We also found that the quark mass difference 
and the electric charge neutrality (but not the color neu- 
trality) play key roles in determining the phase struc- 
ture near the critical temperature. In weak coupling, 
it was shown that the following successive phase tran- 
sitions occur near the super-normal phase boundary: a 
modified color-flavor locked (mCFL) phase (ud, ds, and 
us pairings) — > a "dSC" phase (ud and ds pairings) — > 
an isoscalar two-flavor superconducting (2SC) phase (ud 
pairing) — > a normal phase (no pairing). 

The purposes of this paper are (i) to give a detailed 
account of the results given in Ref. [HI and (ii) to in- 
vestigate the elementary excitation modes (gluons with 
the Meissner mass, massless gluons, gapped quarks, and 
gapless quarks). 

The content of this paper is as follows. In Sec. II, we 
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start with a toy model that captures the essential features 
of the hierarchical structure of the phase transitions. 
Then we construct a general form of the GL potential 
with the quark masses in such a way that it fulfills sym- 
metry constraints. In Sec. Ill, we evaluate the param- 
eters characterizing the GL potential terms in the weak 
coupling region with unequal quark masses [mi, i=u, d, s) 
and unequal quark chemical potentials (//», i=u,d,s) by 
utilizing the Cornwall- Jackiw-Tomboulis (CJT) effective 
action [Ijj. In Sec. IV, we consider a simplified situation 
( m u.d — and m s ^ 0) and calculate the corrections to 
the GL potential from nonzero m s . In Sec. V, we intro- 
duce a parametrization of the gap and analyze possible 
phases near the super-normal phase boundary at finite 
temperature. For these phases, residual symmetries, the 
Meissner masses of gluons, and the gapped and gapless 
quark modes are also investigated. Section VI is devoted 
to summary and concluding remarks. 



II. GINZBURG-LANDAU POTENTIAL 

In this section, as an introduction to later sections, we 
first study a toy model that captures the essential fea- 
tures of such a hierarchical structure of phase transitions 
as encountered in dense quark matter. Validity of the 
various approximations adopted later will be also clari- 
fied by the analysis of the toy model. We then construct 
a general form of the GL potential with the quark masses 
such that it satisfies symmetry constraints. 



we parametrize these coefficients as 
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where \i and T c are the quark chemical potential and the 
critical temperature in the chiral limit, m and g are the 
quark mass and the QCD coupling constant, respectively. 
In the last equalities of Eqs. we mimic the ac- 

tual parametric forms of t, a, and shown in Eqs. i|28|) 
and H61fl except for numerical factors. We do not have to 
consider the quark mass corrections to fa and fa under 
the condition to be shown shortly. 6 is a key parameter 
which governs the multiple structure of the phase transi- 
tions: Due to the presence of 5, the temperatures where 
the coefficients affixed to X 2 and Y 2 change sign are no 
longer identical. 

The minimum of the potential can be simply obtained 
from dVL/dX = 80,/dY — 0. Then one finds three differ- 
ent phases: 

1. a phase where both X and Y have non- vanishing 
condensates, 

X 2 = -±( a + 38) = -%(t + 3v), (5) 



Y 2 



1 

60 



(a - 35) 
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(t-3a). 



(6) 



A. A toy model 

Let us assume two real order parameters, X and Y, 
and consider a GL potential of the following form: 



2. a phase where only Y has non- vanishing conden- 
sate, 



X A 



Y z 



0, 



1 T 2 

(a-S) = -(t-a). 

Afa ' 4 v ' 



(7) 
(8) 



fi = a(X 2 + Y 2 )+S(X 2 -Y 2 ) 

+fa(X 2 + Y 2 ) 2 + fa(X 4 + Y 4 ). (I) 

This is the most general form of the potential up to the 
quartic order with the "parity" symmetry, X «-> —X and 
Y <-> —Y. For the stability of the potential, fa + fa > 
and 2/3i + fa > are implicitly assumed. 

As we will see later, the coefficient a corresponds to 
the reduced temperature in the chiral limit, while 5 de- 
notes the degree of flavor symmetry breaking (in the 
present context, the breaking of X 2 Y 2 symmetry) 
which is proportional to the quark mass squared. The 
coefficients fa and fa are the quartic couplings that do 
not change sign around the critical point, but more or 
less receive corrections from the quark mass. Bearing in 
mind the results that will be obtained in later sections, 



3. a phase where there are no condensates, 

X 2 = Y 2 = 0. (9) 

The behavior of the condensates as a function of the re- 
duced temperature t is illustrated in Fig. ^ The multiple 
structure of the phase transitions induced by a shown in 
this figure is a general feature which survives in a realistic 
situation that will be discussed in later sections. 

The questions to be addressed here are (i) what is the 
region of temperature where the GL potential Eq. (JTJ 
expanded up to the quartic order is valid, and (ii) what is 
the justification of neglecting the quark mass corrections 
to fa and fa. These questions can be answered if we 
restrict ourselves to the region 

|i| < const • a. (10) 
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FIG. 1: Condensates X 2 and Y 2 as a function of the reduced 
temperature t in a toy model, a is the parameter that breaks 
"flavor" symmetry X 2 — > Y 2 . 

Within this interval, both X 2 and Y 2 are of order or even 
smaller than oT 2 . Therefore, the dimensionless expan- 
sion parameters of the GL potential satisfy 



(X 2 ,Y 2 



<0(a), 



(11) 



which in turn guarantees that the higher order terms be- 
yond the quartic order are negligible. Furthermore, the 
additional quark mass corrections to (3\ and /3 2 give only 
higher order corrections to X 2 and Y 2 as can be seen 
from, e.g., Eqs. (0 and JfjJ); X 2 and Y 2 are already pro- 
portional to a if Eq. (|I0(l is satisfied. 



B. General form of GL potential 

In this subsection, we construct a general form of the 
GL potential on the basis of the QCD symmetry under 
G = SU(3) C x SU(3) L x SU{3) R x U{l) B - 

For this purpose, let us first consider a basic variable, 
the gap function 4>bcjk- It is defined through the pairing 
gap A for a quark (color 6 and flavor j) and another 
quark (color c and flavor fc) on the Fermi surface as will 
be shown in Sec. Ill [in particular, see Eqs. (|18|1 and 1)30(1 ], 
Focusing our attention on spin-zero pairings of positive 
energy quarks in anti-symmetric combinations of colors 
and of flavors, we may rewrite 4>bcjk as 



\L,R- 



(12) 



Here the 3x 3 matrix [c^Jl belongs to the fundamental 
representation of SU(3)c x SU(3)l- Similar properties 
hold for [d l a ]ji under the change R <-> L. In this subsec- 
tion, we omit the color and flavor indices and write the 
gap function as dz,^ for simplicity. 

Under the operation of Gl.r = SU(3)c x SU(3)l,r x 
U(1)b, the left- and right-handed quarks and the gap 
di p transform as 



&L,R 



e lv U c U L MiP L ,R, 
e- 2l *U L , R &L, R Uc, 



(13) 
(14) 



where Ul,r?s are 3x3 unitary matrices corresponding to 
each SU(3)l.r symmetry. The phase factor is associated 
with the U(1)b rotation. 



We now try to incorporate small but nonzero quark 
masses and write down possible terms of the GL potential 
allowed by the QCD symmetry G. In QCD with quark 
masses, we have such a term as 



C = ipLmipR + h.c. 



(15) 



with m being a 3 x 3 mass matrix in flavor space. Just 
like the standard procedure to construct the chiral La- 
grangians, m is assumed to transform in the following 
way Q, 



m — > U LiTiUl 



(16) 



Following the transformation laws, Eqs. (|14(l and i|16[l . 
we can write down a possible form of the GL potential 
as an expansion in terms of the order parameter, which 
is, in the present context, the on-shell gap function d^ n. 
Following the discussion in Sec. II. A, we take into account 
the quark mass terms up to 0(m 2 ) only in the quadratic 
terms in dL,R and obtain 

n = (a ) L Tr(dld L ) + (ao)flTr(dJjdfl) (17) 
+ai [Tr(d^mdfl) + h.c] 

+(ffl2)LTr(d[mm , d L ) + (a 2 )flTr(d^m t md 7? ,) 

+(a2)LTr(d^di)Tr(m t m) + (a 2 )_RTr(d^,dfl)Tr(m t TO) 

+c[detm • Tr(d^ R m~ 1 d L ) + h.c] 

+ (6 1 ) i [Tr(d[d i )] 2 + (6 1 ) fl [Tr(dt ? d fl )] 2 

+ (6 2 ) i Tr[(dld i ) 2 ] + (6 2 )«Tr[(d^) 2 ] 

+6 3 Tr[(d[d i )(dl ? dfl)] + 6 4 Tr(d[d L )Tr(dt ? d«), 

where the a's, 6's, and c are the expansion coefficients 
being functions of temperature and chemical potential. 

Only the term proportional to a\ in Eq. (|17f) breaks 
U(1)a symmetry and thus originates from anomaly in- 
duced interactions. It does not affect the phase structure 
near the super-normal phase boundary at asymptotically 
high density as will be shown in Sec. IV. D. The term pro- 
portional to c also does not play an important role near 
the phase boundary as will be discussed at the end of 
Sec. III.B. 

We can further take into account the effects of electric 
and color neutralities and discuss possible terms to be 
included in Eq. (|17fl . In particular, the electric charge 
neutrality plays a central role in this paper. Instead of 
writing down its general structure, we will show its ex- 
plicit form calculated in the weak coupling regime in Sec. 

III. The effect enters into the quadratic terms of the GL 
potential. As for the color neutrality, its effect is sub- 
dominant in the vicinity of the phase boundary as has 
been already shown in Ref. @; the color neutrality af- 
fects only the quartic part of the GL potential (see Sec. 

IV. C). 
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III. CORRECTIONS TO THE GL POTENTIAL 
IN WEAK COUPLING 

In this section, we calculate corrections to the GL po- 
tential by unequal quark masses (mj, i=u,d,s) and un- 
equal quark chemical potentials {pi, i=u,d,s). To sim- 
plify the derivation, we make several ansatze for the pair- 
ing gap: 

1. The pairing between quarks is assumed to be in 
the zero total angular momentum (J = 0) chan- 
nel 0|- Furthermore, it is assumed to be in the 
LL and RR channel as in the massless case and 
in the positive parity channel as in the presence of 
U(1) A breaking @. Then the gap can be written 
as A(k) = 7 5 AW + 7 . k 7 ° 7 5 A^ Q. 

2. The pairing is projected onto positive energy states. 
In this case, it is convenient to rewrite A(fc) as 



A(fc)=7 5 ^ ,k)A+(k) 



(18) 



where A + (k) = (1 + 77 • k)/2 is a projection op- 
erator onto the positive energy states of massless 
quarks. 

3. The pairing is assumed to take place in the color 
antisymmetric and flavor antisymmetric channel. 
This color channel is indeed the most attractive 
in weak coupling 0, 0, El , while the flavor chan- 
nel is chosen in such a way as to satisfy the Pauli 
principle. Then, the pairing gap of a quark (color b 
and flavor j) and another quark (color c and flavor 
k) at the Fermi surface takes the form as given in 
Eq. (|12[) . Furthermore, the pairing in the positive 
parity channel is represented by 

d L =d R = d. (19) 

4. We adopt the fact that the three-momentum de- 
pendence of the on-shell gap A(q = e(q),q) or, 
equivalently, 0(q;T) = <j){ck\ = e(q),q) can be ap- 
proximately factorized as 01 



(q;T)~0(|q|= W r)/(q), 



where 

/(q) 
9 

b 

y(q) 
£(q) 



= In 



.9/3V27T, 
256^ 4 (2/3 5 2 ) 5 / 2 , 
26/i 



|q|-/i|+£(q) 



/■'■■I 



I0(q;r)| 



21 1/2 



(20) 

(21) 
(22) 
(23) 

(24) 
(25) 



5. We will assume later in Sec. V that d l a is diago- 
nal in color-flavor space even under the influence 
of nonzero quark masses and charge chemical po- 
tentials. Then we will minimize the GL potential 
within the subspace. 



In the following, we start with a known result for mass- 
less quarks (Sec. III. A) and then proceed to describe a 
more general situation with unequal quark masses (m*, 
i=u,d,s) and unequal quark chemical potentials 
i=u,d,s) in Sec. III.B. 



A. Case with massless three flavors 

For a homogeneous system composed of massless 
quarks (m u ^d, s = 0), quarks have a common chemical 
potential /x, and the GL potential near the critical tem- 
perature T c expanded up to quartic order reads 0, Il8| 

n = 5^ |d a | 2 + Kl 2 ) 2 + E l d °* • d " I 2 - 



a b 



(26) 



where (d a ) 1 = (d", d^,d^), and the inner product is taken 
for flavor indices. Using Eq. . one can relate the 
coefficients in Eq. (|17H to those in Eq. (|26[1 as 



Pi = (&i)z, + (bi) R + 64, 
h = {b 2 ) L + (b 2 )n + h, 
ol = (& )l + (ao)/i. 



(27) 



This potential is manifestly invariant under SU (3)c x 
SU(3)l+r x U(1)b rotation. In the weak coupling ap- 
proximation where the one-gluon exchange force in the 
normal medium is responsible for the pairing, the coeffi- 
cients have been calculated as 61 



0i = A 



7C(3) 
8(nT c y 



■N(ji) =13, a = AN(fi)t = a t. (28) 



Here N(fi) = fi 2 /2ir 2 is the density of states at the Fermi 
surface, and t = (T — T c )/T c is the reduced tempera- 
ture. With the parameters in Eq. (|28|l . one finds a single 
second-order phase transition at T = T c from the CFL 
phase (d l a oc S l a ) to the normal phase (d l a — 0) in the 
mean- field theory [(|. 



B. Case with unequal quark masses and chemical 
potentials 

For a homogeneous system composed of unequal quark 
masses, there arise differences in the chemical potential 
among flavors and among colors due to charge and color 
neutrality conditions. The GL potential in this case is 
obtained in a similar approach to that adopted in Ref. 
|3| where a correction to Eq. (|26l) from color neutrality 
is calculated. 

First we start with the Nambu-Gor'kov two compo- 
nent field (ipaii^ai)' wnere ^° = CV> T is the charge- 
conjugate spinor. The quark propagator of this field with 
the Hartree-Fock contributions ignored in the diagonal 
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part reads 



G(k) 



( G {11 \k) G (12) (fc) 
\G {21 \k) G^ 2 \k) 

jk + j°A4 — m 

A(k) jk 



A(k) 
1 °M T 



(29) 



(30) 



Here A = 7°A t 7°, and 



M-abij — <>ab<>ijHi, m a bij = Sab^ijlTli, 



(31) 



are the quark chemical potentials and the quark masses in 
color (a, b) and flavor space, respectively. As in Sec. 
III. A, we again define ^ as a quark chemical potential in 
the chiral limit (m Ut d. s — 0). Once unequal quark masses 
are included under charge neutrality and (3 equilibrium, 
Hi differs from h- We define the deviation of A4 a uj from 
the massless case as 5M. a bij = 8 ab 5ij{Hi — A*)- 

The free quark propagator Go and the self-energy X are 
defined from the diagonal and off-diagonal components of 
G -1 as usual: 



G _1 (fc) = G ( 7 1 (fe)-£(fe) 



(32) 



Then the GL potential, which is a difference between 
the superfluid and normal phases near the critical tem- 
perature, is written in the CJT form |T^ 



tt = n 1 + n 2 

f d 3 q 
2 (2tt)3 



(33) 

''" 9 Tr[-G(q)E( q )+lnGo 1 (q)G(q)}, 



O - ^T 2 ^ f d 3 k f d 3 q 



(34) 



x Tr 



VffiQ ~ k)^-G^(k)^ - G^(q) 



+V^{ q -k)^[ — \ G^(k)^-G^(q) . (35) 



Here VL2 represents the two particle irreducible graphs in 
the mean-field approximation with T>(q) being the gluon 
propagator in hard dense loop approximation, and A" the 
generators of color SU(3). The summations are taken 
over Matsubara frequencies of quarks. 

In the following, we will consider the corrections from 
S Mabij and ra a uj to the quadratic term in A assuming 
that the corrections are small. The corrections to the 
quartic terms are negligible near the critical temperature. 
Figure 2 represents the C(A 2 ) contributions to f^i and 

n 2 . 



First we consider 



n odd 



1 s 

d 3 



Tr 



1 



79 + 7°A1 — m 
A(?) 



A(<?) 



El 



abij 



79 — 7°A1 T — m 

- — %-r tanh [ -^= \ Tr 
(2tt) 4 \2TJ 



jq + 7°/ii - mi 



x A(q) abij 



1 



79 - 7 ft 



-A(g)i 



.(36) 



The first quark propagator in the right-hand side of Eq. 
(|36[) has particle and anti-particle poles, 

eP ± = ±H~p F , (37) 

while the second propagator has hole and anti-hole poles, 

(38) 



e hT = T|q| 



■Pf, 



where p F — y Hi — mf is the Fermi momentum of flavor 
i. We define the shift of the chemical potential of flavor 
i from /i as 

Sfj,i = Hi - /i. (39) 

Then, up to 0(m 2 / fx 2 , Sfi/ fJ-), Pp can be written as 

p F ~ m - m 2 /2fi. (40) 

Of the poles in Eqs. (|37|) and l|38|l . only two of them, 
q a = e p+ (particles) and q = e h ~ (particle holes), have 
relevant contributions to the potential. The other two 
poles contribute to the gaps on antiparticle (antiparticle 
hole) mass shell, which are irrelevant in our calculation. 
Let us define 



q = |q| - M> 

tanh(g/2T) 



m = 



9 



(41) 
(42) 



Then, combining the residues of the two poles, we obtain 
d 3 q 



Here. 



E 1 
^ 2 

abij 



r2/T]l v F( 5 )-^^)|„,,„ y |^. N :j, 



Spij = (Pf+P^)/2-M- 



(44) 



For later purpose, we also introduce an averaged shift of 
the chemical potential as 



(45) 



In the diagrammatic language, the first term in the 
right-hand side of Eq. (|43|l corresponds to (aO) in Fig. 
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|5J The second term corresponds to (al) and (a2). (a3) 
vanishes because of the Lorentz structure of the gaps to- 
gether with the identity, Tr(A+A~) = 0. 

To perform the momentum integration in Eq. (|43[) . it is 
convenient to introduce a positive constant k ~ 0(<7 _1 ) 
which cancels at the end of the calculation . The 
cutoff of the three-momentum A is taken to be bfj, for the 
sake of consistency with the formula, Eq. I|21|) . We divide 
the region of the momentum integral —fjt < |q| — /J, < 
b\x into two parts, < ||q| — fi\ < k4>q (region I), and 
— fi < |q| — fx < —k4>o and k</>o < |q| — fJ, < bfi (region 
II) [13. Here we defined <f> = [(1/9) ^ \4>aUj {n, T = 

abij 

0)| 2 ] 1/2 - In each region, the following approximations 
can be safely taken: in region I, smgy\T=o ~ L and in 
region II, tanh [(|q| - M )/2T] ~ (|q| - M )/||q| - 

The first term in the right-hand side of Eq. 143(1 is 
calculated as 



(aO) 



abij J y ' 

1 { r K ^° 

abij 1 I J( > 



-lnb 



\ d yf 2 (y) \\^abij\ 



(46) 



In the first equality, we replaced the momentum in the 
measure by the density of states at the Fermi surface. In 
the second equality, we used the weak coupling approx- 
imation g <C 1 and expanded f(y) in terms of g. We 
also used the explicit g dependence of T c and 4>o in weak 

coupling [SEI, !n(T c //x) - ln(0 o /^) 3x 2 /V2g, to 

obtain the final form. 

The second term in the right-hand side of Eq. (|43|) is 
more involved: 



(al) 



(a2) 
1 

~2 

1 

4^2 



d 3 q dT{q) 



dq 



6 P M ! 



(2tt) 3 

bfj, 

dq {q + ^i) 2 {T{q)-T{q 



where we abbreviated internal indices for simplicity. 
f(bfi) appearing in the first term of the last bracket in 
Eq. is 0(g(4> 2 /(b/i) 2 )) and can be neglected at high 
densities. Then the second term leads to 



Sp)}f(q) 



dqT{q)f\q){q + ») 




dq T{q)f 2 (q){q + ^Y 



(al) + (a2) 




dq f(q) 



phi b 




/ + 


f )dyf 2 (y)\ 


/ \n(bfi j K,<fio) 


Jln(bp,/ K <p ) ) J 



(48) 



Here we again expanded f(y) in terms of g and used the 
explicit g dependence of T c in weak coupling to obtain 
the final form. 

Then the sum of Eqs. I|46|l and 1|48[) becomes 



ftf = (aO) + (al) + (a2) 



X>(m) 

abij 



In 



5p 



13 In ( ^ 



T 

\4>abij\ 2 ■ 



(49) 



A 2 

Let us next evaluate VL 2 ■ In the diagrammatic de- 
scription, this corresponds to (b0)-(b4) in Fig.[5]in which 
the quark propagator is expanded in powers of masses 
and chemical potential differences. Of these graphs, (b2) 
vanishes because of the Lorentz structure of the gaps, 
and (b3) vanishes because the gaps are constructed by 
quarks having the same chirality. As a result, only (bO), 




Sp h 5 



FIG. 2: The quadratic terms of the GL potential, x and + 

A 2 

represent m and Sfi, respectively. (a0)-(a3) represent Q 1 , 

A 2 

and (b0)-(b4) are the two particle irreducible graphs Q 2 ■ 
(aO) and (bO) are the quadratic terms of the potential in the 
m = Sfi — case. Among these graphs, (a3) and (b2) vanish 
due to the Lorentz structure of the gap. (b3) vanishes because 
we consider the gaps constructed by quarks having the same 

/ chirality. Only (al), (a2), (bl), and (b4) give finite correction 

^ ' to the original GL potential [(aO) plus (bO)]. 
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(bl), and (b4) have a finite value. We expand these dia- 
grams by t = \T - T c \/T c around T c up to 0(t) for (bO) 
and 0(1) for (bl) and (b4). 

We use a relation between T c and g when rrii and 5pi 
in the denominator of the quark propagator are taken to 
be zero 0, 



/(q) 



9 2 J (2ttH 



d 4 k , / fen , 

tanh — Tr 
2T r 



x ( y) T ^ " V^)"^* + 7°A*r V ( y ) ./ (ki 



(50) 



This equation is derived by taking T — > T c in the gap 
equation with the decomposition Eq. 1)20(1 . 

Using Eq. (|50|l which relates T c and g, we can evaluate 
(bO), (bl), and (b4) in Fig.0 We derive (bO) as follows. 
In leading order in t, (bO) becomes 



(bO) 



xTr 



d 4 q , , / go 

- — — tann 

2tt 4 I 2T r 



0(1(1 



^(»-%"(t 



xG< 2 »»)7''(y))g ,12) (9) 



T=T C 



^. tanh f — 1 
(2tt) 4 V2T/ 



x Tr 



/(q)0(T)G( 12 )(g)' 



T=T C 



^ tMi (^) Tr[ A(,) G <")(,)l 



T=T C 



-fi^ |t=t c 



The 0(t) term of (bO) becomes 
(bO)- 

(T - Tc) 4dr{/(2^ 



(51) 



2T 



tanh £ Tr[G( 21 )(fc)A(fc)] 



+ 



T=T C 



2(T-T C ) — (-fif )| r=Te 



(52) 



(bl) and (b4) can be calculated in a similar way, and 
they become 



(bl) + (b4)=7V( M )^ln(^) \<f> ab 



(53) 



Adding (bl) and (b4) to (bO), we finally obtain 
flf = (bO) + (bl) + (b4) 



= -nf\ T = Tc + 2(T-T c )^(-nf)\ T=Tc 

+ N(fi)^-ln(^) \<f> abij \ 2 



X \4>abij\ 



(54) 



Combining Eqs. (|49|l and 154() . we obtain the final re- 
sult for the quadratic term of the GL free energy, 

■«-£MK£)-*H£)} 

X \(j)aUj\ 2 , (55) 

<5p y - ~ S^j - (m 2 + m 2 )/(4fi). (56) 
It can be rewritten in a clearer form, 



GL 

rpij 
_ c 



E 

abij 



= l 



AT( M ) 



T - TV 
T r 



\(/>abij\ 



(57) 
(58) 



Several comments arc in order here. 



• The free energy correction depends only on the av- 
erage shift of the Fermi momenta of paired quarks, 
Spij, as can be seen from Eq. I|55|) . Furthermore, 
the correction may be decomposed into a sum of the 
effects from the averaged chemical potential and 
the averaged mass squared when Sp^ is small, as 
shown in Eq. I|56|) . 



• Equations 157|) and l|58|) indicate that the shift Sp^j 
affects the free energy correction through the crit- 
ical temperature T* J for the pairing. The 
larger the shift Spij is, the higher the melting tem- 
perature Tl 3 becomes. In other words, one can 
compare the thermal stability of two different pairs 
only by taking note of a difference in the average 
density of states between them. 

• The unlocking at T = is different in mechanism 
from our unlocking near T c . At T — 0, color- flavor 
unlocking due to the Fermi momentum mismatch 
between paired quarks is expected at [i ~ m 2 /</>o 
fl9| . In contrast, the Fermi momentum average, is 
important for our unlocking near T c . 

• We can see that in the correction obtained above 
quark flavors do not mix with each other. As long 
as we consider the pairing of positive energy quarks 
and neglect that of antiquarks, the flavor structure 
of a gap always enters in the form of m 2 instead of 
m l m j (i ^ j). 
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IV. GL POTENTIAL IN REALISTIC QUARK 
MATTER 



Up to now we constructed the GL potential for general 
quark masses and chemical potential shifts assuming that 
the corrections are small. In this section, we will focus 
on the case close to a realistic situation by taking m s =/= 
with m u ,d = and by imposing [3 equilibrium and charge 
neutrality conditions. This analysis clarifies the role of 
the strange quark mass in possible phase structures near 
the super-normal phase boundary at finite T . 



Direct m s correction through the quark mass 
term 



Let us first consider the correction directly propor- 
tional to m 2 in Eqs. (|55|l and (|56|l . It is easy to see from 
the flavor structure of Eq. (|55J) that m s affects only us 
and ds pairings. This is also reasonable from the physical 
point of view as will be shown below. 

By inserting m s ^ with m Ui d — and using Eq. 
(|T2|l with Eq. ijTTHl . the direct m s correction in Eq. 
becomes 



Here the coefficient e reads 
m 2 s 

e " a °v ln 



2a Q a, 



(59) 



(60) 



where = 4iV(/i) is defined in Eq. (|28|) . and a dimen- 
sionless parameter a is introduced as 



V8V2/ gv 2 



(61) 



In Eq. HfiOjl we have used a weak coupling relation, 
ln(T c //i) ~ — 37r 2 / (V2.9), which originates from the long- 
range color magnetic interaction [iH Il7l | . 

Since the finite m s decreases the Fermi momentum of 
s quarks, 5p us and 5pd s are smaller than Sp u d, and thus e 
becomes positive such that ud pairing is favored over us 
and ds pairings. Consequently, the CFL phase becomes 
asymmetric in flavor space and its critical temperature 
is lowered, leading to the appearance of the 2SC phase 
(4 cx 5 ts ) just below T c H|. 



B. Indirect m s correction through the charge 
chemical potentials 

We proceed to discuss the correction from the charge 
chemical potentials, which is proportional to Sfiij in Eqs. 
(|55|l and (|56|) . Under /3 equilibrium and charge neutrality, 
the electron chemical potential /i e and the shift 6 fit of the 
chemical potential of flavor i are related as 



with qi being the electric charges of the quarks; q u = 2/3 
and qd, s = —1/3. In weak coupling, where one may re- 
gard normal quark matter and electrons as noninteract- 
ing Fermi gases, /i e is related to m s as 



(63) 



This estimate is valid in the vicinity of T c where correc- 
tions to fi e by a finite pairing gap affect only the quartic 
terms in the GL potential. 

By inserting Eqs. (|62|) and (|63|l and using Eq. with 
Eq. (|I9() . the indirect m s correction in Eq. I|55(l becomes 



^Ei d °i 2 -EKi 2 )- 



Here the coefficient 77 reads 

m 2 s ( fi 
^ ao 8^ ln U 



(64) 



(65) 



Since 5 /ids is larger than 8[i u d and Sfi us and it increases 
the average Fermi momentum of ds quarks in Eq. l|56|l . 
the indirect m s correction favors ds pairing than ud and 
us pairings. 

Note that we only consider the modification of T c due 
to nonzero m s through the properties of the superfluid 
phase. Actually, T c is modified not only by the properties 
of the superfluid phase but also of the normal phase. 
The modification by the normal phase enters through 
the Debye mass in the gluon propagator T> in the relation 
between T c and g, Eq. (|50|l . Fortunately the modification 
due to nonzero m s and charge neutrality in normal quark 
matter is like T c — > T c (l + O(ga)). This modification 
is of higher order than that in superfluid quark matter, 
T c T c (l + 0(a)), to be derived in Sec. V. 



C. Effect of color neutrality 

We consider color neutrality of the system as well. In 
contrast to the case at T = 0, however, it affects only the 
quartic terms in the GL potential because the possible 
chemical potential differences between colors are j(| 0] 



In 

9/x 



7 )! 



(66) 



(62) 



which are already of quadratic order in the gap. Here 
S^ab = [tab — H&ab, and [x a b is the chemical potential con- 
jugate to n ab = E ^6i7°^ai 111- In weak coupling the 

i 

magnitude of the correction to the quartic terms is sup- 
pressed by 0((T C / gy) 2 ) compared to the leading quartic 
terms. Thus color neutrality has no essential consequence 
to the phase transitions considered in this paper. A ma- 
jor difference between the corrections from the charge 
neutrality and the color neutrality is that the former 
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shifts the quark chemical potentials even in the normal 
phase, while the latter works only when the pairing oc- 
curs. This is why the former is more important than the 
latter near T c . 



D. Effect of instantons 

In QCD with three flavors, instantons induce chirality- 
flipping six-fermion interactions. This interaction leads 
to a sextic term in the gap in the chiral limit, which 
is irrelevant near the super-normal phase boundary. If 
the strange quark mass enters, the direct instantons in- 
duce an effective four-fermion interaction between u and 
d quarks [22|. This leads to a quadratic term in the GL 
potential, £ Ma 1 2 > which corresponds to the a\ term in 

a 

Eq. H17JI . An explicit calculation in weak coupling shows 
that £ ~ -a {m s /ii){kQ CI >liif(\/g) 14 . The negative 
sign indicates that the instanton effect favors ud pair- 
ing as does one-gluon exchange [see Eq. i|59fl]. However, 
the magnitude of £ is highly suppressed at high densities. 
Therefore we will ignore this term in the following. 

In summary of Sees. IV.A-D, as far as we are close 
to the super-normal phase boundary at high density, we 
have only to consider the corrections from the strange 
quark mass and the electric charge neutrality, which favor 
ud pairing and ds pairing, respectively. 



E. Validity of the approximations 

Our analysis of the GL potential near T c is valid as 
long as (a) g < 1, (b) a < 1, and (c) T c < gfi. The 
condition (a) is necessary for the dominance of the long- 
range magnetic interaction acting between quarks. The 
condition (b) allowed us to consider the corrections to the 
GL potential by the strange quark mass and differences 
in the chemical potential among flavors up to quadratic 
order in the gap. The condition (c) is imposed so that the 
correction to the quartic terms in the GL potential from 
the color neutrality can be neglected. These conditions 
are all satisfied at asymptotically high density. 

If we use the result for T c calculated in weak coupling 
EE 113, one finds that (c) is a consequence of (a). The 
conditions (a) and (b) can also be combined into 

o 

mz 

-f « g « I. (67) 



V. MELTING PATTERN OF DIQUARK 
CONDENSATES 

In this section, we clarify the phase structure near the 
transition temperature using the GL potential corrected 
by Eqs. I|59|l and l|64|) . We start from assuming the form 



of the condensation and derive the temperature depen- 
dence of the gaps that minimize the potential. Then we 
discuss the properties of quark quasiparticles and trans- 
verse gluons in each of the phases near the transition 
temperature. 

A. Phase structure with finite m s and charge 
neutrality 

Since the two effects of nonzero m Sl characterized by 
Eqs. (|59|l and l|64|) . favor ud pairing and ds pairing, 
respectively, the finite temperature transition from the 
CFL to the normal phase at m s = is significantly mod- 
ified. In fact, as we show in detail below, successive color- 
flavor unlockings take place instead of a simultaneous un- 
locking of all color-flavor combinations. To describe this 
hierarchical thermal unlocking, it is convenient to intro- 
duce a parameterization 

/Ai \ 
4=0 A 2 . (68) 
V A 3 J 

We assume ^1,2,3 to be real. We also name the phases 
for later convenience as E3 

Ai, 2 .3^0 : mCFL, 

Ai - 0, A 2 , 3 + : uSC, 

A 2 = 0, Ai,3^0 : dSC, (69) 

A 3 = 0, Ai 2 ^ : sSC, 

Ai, a =0, A 3 ^0 : 2SC, 

where dSC (uSC, sSC) stands for superconductivity in 
which for d (it, s) quarks all three colors are involved in 
the pairing. 

In terms of the parameterization (|68|l . the GL potential 
with corrections of 0(m 2 s ) to the quadratic term, Eqs. 
IJ55JI and (H2J, reads 

n = at' (&t + A 2 . + A3) - eAg - 77A 2 

+ ft(A? + A 2 2 + A3) 2 + (3 2 (A{ + A 4 2 + Aj), (70) 

where a' = a + £ H — ■ 
3 

We proceed to analyze the phase structure dictated by 
Eq. (|70l with the weak coupling parameters |[BSJ, 
and H65|) up to leading order in g. A stable condensation 
must satisfy the gap equations 

dQ/dAt <2 , 3 = 0, (71) 

and it must minimize the potential Eq. I|7UI) . We compare 
the energies of all condensates in Eq. H69[) and decide the 
most stable condensation at each temperature near T c . 

In Figs. 3 and 4 the results thus obtained for the phase 
structure near T c are summarized. 

Figure 3(a) shows the second-order phase transition, 
CFL — > normal for m s — 0. Figures 3(b,c) represent how 
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CmCFL> 

(a) (b) (c) 

FIG. 3: Transition temperatures of the three-flavor color su- 
perconductor in weak coupling: (a) all quarks are massless; 
(b) nonzero m s in the quark propagator is considered; (c) 
electric charge neutrality is further imposed. The numbers 
attached to the arrows are in units of aT c . From Ref. 



ered, and the 2SC phase is allowed to appear at temper- 
atures between T* and T c . As one goes from (b) to (c), 
the average chemical potential of ds (ud and us) quarks 
increases (decreases). Accordingly, the transition tem- 
peratures further change from T c to T* 11 and from Xjf to 
T] and T c n . 

Now we examine in more detail how the color-flavor 
unlocking in case (c) proceeds with increasing T from 
the region below Tl 

(i) Just below r], we have a CFL-like phase, but the three 
gaps take different values, with an order A3 > Ai > 
A2 7^ (the mCFL phase). The reason why this order is 
realized can be understood from the GL potential 17UII . 
The e-term and ry-term in Eq. (|70|l tend to destabilize us 
pairing (A2) relative to ud pairing (A3) and ds pairing 
(Ai), and since e > r](> 0), ds pairing is destabilized 
more effectively than ud pairing. The value of each gap 
in the mCFL phase reads 



mCFL dSC 2SC 
► < 




1c 1c 



A^ 



op 

8/3 



T r -T 



+ Ci(J 



(72) 



with ci j2: 3 = (-4/3,-16/3,8/3). The mCFL phase has 
only a global symmetry U (1)c+l+r X U(\) c +l+r in con- 
trast to the global symmetry SU(3)c+l+r in the CFL 
phase with m Ui d, s = 0. The remaining U(l) generators 
are 



Qi 

Q,2 



T 3 



S3, 
S$. 



(73) 



FIG. 4: A schematic illustration of the gaps squared as a 
function of T. From Ref. [TH|. 



the phase transitions and their critical temperatures bi- 
furcate as we introduce (b) effects of a nonzero m s in the 
quark propagator and then (c) effects of charge neutral- 
ity. In case (b), two second-order phase transitions arise: 
mCFL (with Ai = A 2 ) -> 2SC at T = T c s = (1 - 4ct)T c , 
and 2SC — > normal at T — T c . In case (c), there arise 
three successive second-order phase transitions, mCFL 
-► dSC at T = T l dSC -» 2SC at T = T c n , and 2SC -> 
normal at T = T™. Shown in Fig. 4 is the T-dependence 
of the gaps ^1,2,3 for the case (c). All the gaps are con- 
tinuous functions of T, but their slopes are discontinuous 
at the critical points, which reflects the second order na- 
ture of the transitions in the mean-field treatment of Eq. 

We may understand the bifurcation of the transition 
temperatures as follows. In the massless case (a), T c is 
degenerate between the CFL and 2SC phases, the chem- 
ical potential is common to all three flavors and colors, 
and the CFL phase is more favorable than the 2SC phase 
below T c . As one goes from (a) to (b), the density of 
states of the s quarks at the Fermi surface is reduced. 
Then the critical temperature for the CFL phase is low- 



Here T a (Si) are the part of the generators of 
SU(3) C (SU(3) L+R ) with explicit forms T 3 = S3 = 
(l/2)diag(l,-l,0) andT 8 = Sg = (l/2v / 3)diag(l, 1, -2). 
The free energy in this phase is 



rnCFL 



3og 
16(3 



T-T c 



3 



(74) 



As T increases, the first unlocking transition, the unlock- 
ing of A 2 (the pairing between Bu and Rs quarks), takes 
place at the critical temperature, 



T!={l-fa)T,. 



(75) 



(ii) For T l c <T < T c n , A 2 = and 



A 2 a {T c -T 



66 



+ CiCF 



(76) 



with ci,3 = (—7/3,2/3). In this phase, we have only 
ud and ds pairings (the dSC phase) , and there is a mani- 
fest symmetry, U(1) c +l+r* U{\) c +l+r^ U(1) c +v+b x 
U(1)c+v+b, where the corresponding U(l) generators 
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= ^3 


Q 

— «J3, 








= t 8 








Q 3 


= Q- 




-25 3! 




Q4 


- Q- 


f 7! 5s 


-25 3j 


(77) 



respectively. Here Q — 2/3 is the generator of baryon 
charge (Q) and acts on the gap as Eq. (|Tlfl . The free 
energy in this phase is 



■dsc = —77-5 




T-T c 
T, 



6 



(78) 



At T = TJ 1 , the second unlocking transition, the unlock- 
ing of Ai (the pairing between Gs and Bd quarks), takes 
place at the critical temperature, 



T 



7 



a T r 



(79) 



(iii) For 7; 11 < T < T c m , one finds the 2SC phase, which 
has only ud pairing with 



(T c -T 1 



Ao = — . 

3 4/3 V T t 



r 



(80) 



The 2SC phase has a symmetry SU(2) C x SU(2) L+R x 
?7(1)c+b x ?7(1)l+_r+b, where the corresponding U(l) 
generators are 



Q 5 = Q + V3T 8 , 
Qe = Q + V3S 8 , 

respectively. The free energy in this phase is 
al (T-T c ( 1 







2SC 



8/3 V T c 



(81) 



(82) 



The final unlocking transition where A3 (the pairing be- 
tween Rd and Gu quarks) vanishes occurs at 



T. 



in 



i--.it; 



(83) 



Above TjF 1 , the system is in the normal phase. Note 
that in the temperature region (TjJ — T)/T* < const • a, 
the gaps satisfy the condition Qllfl. so the GL potential 
expanded up to quartic order in the gap is valid. 

So far, we have derived the parameters /?i,2, e, and 77 
in the weak coupling approximation and discussed the 
phase structure expected at asymptotically high density. 
Now we relax the weak coupling constraint and study 
possible phase structures expected from the GL potential 
of the form Eq. (|70|) with arbitrary coupling strengths. 
To simplify the argument, we assume /3i = 02 > and 



A 3 < A 2 < A1 



sSC type 



A3 < A, < Ai/ 



A 2 <A 3 <Ai , 



dSC type 



/ A 2 <A, <A 3 



uSC type 



7]-e/ 

/ Ai < A 3 < A 2 



Ai < A 2 < As 



FIG. 5: The phase structure obtained from Eq. I70H in (e, rj) 
plane for f3i = P2 > 0. The mCFL phase is assumed to 
appear at sufficiently low temperature. ® indicates the point 
calculated in the weak coupling approximation. 



draw the phase diagram in the space of the parameters e 
and 77 in Fig. [S] 

If we consider the case where the mCFL phase appears 
at temperature sufficiently below T c , the relative magni- 
tudes of the gaps behave as shown in Fig. 03 If A2 is the 
smallest of the three gaps, we call the melting pattern 
"dSC-type" because the dSC phase would appear as T 
increases. "uSC-type" and "sSC-type" are defined in a 
similar way. As in Fig. depending on the relative mag- 
nitude of e and r], the dSC-type melting pattern is divided 
into two different classes in which the 2SC-type phase 
(Ai ^ 0, A 2 , 3 = or A 3 ^ 0, Ai )2 = 0) appears just 
below the transition temperature to the normal phase. 
The uSC-type and sSC-type melting patterns also have 
such a two-fold structure. The hierarchical thermal un- 
locking (mCFL — > dSC — * 2SC — > normal) is realized in 
weak coupling as shown by ® in the figure. The figure 
indicates that it is a rather robust phenomenon near the 
super-normal phase boundary. In other words, the uSC 
or sSC phase appears only when e or r) changes sign at 
least under the form of the GL potential in Eq. (|7C)p with 
an assumption /3i = 02 > 0. We mention here that a re- 
cent analysis in Ref. [23] using the Nambu-Jona-Lasinio 
(NJL) model shows that the signs and the ratios of the 
coupling strengths in the GL potential take similar values 
with our weak coupling values: 0i/02 = 1 (with 0\ > 
and 02 > 0) and C = e/v = 2-4 (with e > and 77 > 0). 



B. Properties of quark and gluon modes 

Next we resume considering about the weak coupling 
limit and discuss the properties of quark and gluon modes 
such as the excitation energies in each of the mCFL, dSC, 
and 2SC phases. The low energy excitations of these 
phases are massive or massless gauge bosons, gapped or 
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gapless quark quasiparticles, and modes that are associ- 
ated with fluctuations in the diquark fields. Here we con- 
centrate on the gauge bosons and gapless fermions. As 
for gauge bosons, we only consider the transverse com- 
ponents of the gauge fields which, in the static limit, can 
become massive only in a superconductor. The longitudi- 
nal components are Debye screened in the vicinity of the 
transition temperature. Note that the "gapless fermions" 
in our previous paper ^lj have been defined only by di- 
agonalizing the self-energy £ in Eq. I|32l) . In the present 
paper, we have diagonalized the full propagator G^ 1 (k) 
to define the gapless fermions in a more precise manner 
as will be shown below. 

In Table I, we summarize the symmetries, the num- 
ber of massive gauge bosons, and the number of gapless 
fcrmion modes in each phase. As we will show in detail 
below, the number of gapless quark modes depends not 
only on which phase we consider but also on the temper- 
ature itself. In each phase, in addition to the unpaired 
quark modes whose excitation energies are naturally zero, 
gapless modes appear if the temperature is in the region 
where the gaps satisfy the conditions listed in Table I. 

In each phase, the gauge symmetry is broken partially 
or even totally. The Meissner masses for the general 
condensate, Eq. and the corresponding massive (or 
massless) gauge fields can be calculated using the GL 
theory coupled to gluon fields as 





= m A 2 = 
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(84) 
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Here A is a mixed field of A3 and A* 
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(85) 



Here we followed the notations in Refs. 0,^2- k t is the 
stiffness parameter which controls the spatial variation 
of the gap. In weak coupling it is proportional to the 
quartic coupling (i as 



k t = 13/3 



7C(3) 
24(ttT c ) : 



-jV(/z). 



(86) 



which the full propagator diverges: 
detG _1 (£', q) = G ra.GuG r s ,BuGgs,b<iG R,u,Gd,Bs = 



(87) 



The gap matrix S(fc) can be written in block-diagonal 
form if we adequately take the base whereas the nonin- 
teracting quark propagator Go(fc) is diagonal. The de- 
terminant Eq. H87f) can be decomposed into three 4x4 
determinants Gr<i,Gu, Gr s ,Bu, and Gcs.Bd, and one 6x6 
determinant Gru,g<i,Bs, which may effectively be written 
as 



GRd,Gu = det 
x det 



Gr s ,Bu = det 
x det 



x det 



GRu,Gd,Bs 



det 
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(89) 



(90) 
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(91) 



Here p l F = pp"fo, and pp are the Fermi momenta of i 
quarks. Only in Eqs. (|HHJ we abbreviate the Lorentz 
structure of the gap for simplicity; here A, indicates 
75A + Ai, and Aj indicates 7o(7sA + Ai)7o. Moreover, we 
include the effect of m s as a shift of the effective chemical 
potential p l F following Ref. 25]. Using Eqs. gjfl) , $2$, 
and l|63|l . pp's can be expressed as 



p F = M -(2/3) Me , 
p d F = n+(l/3)n e , 
Pf = M - (5/ 3 )Me- 



(92) 



We can reproduce the results in Table I in Ref. |24| by 
taking a limit Ai,2 — > for the isoscalar 2SC phase and 
Ai = A 2 = A 3 for the CFL phase. 

As given in a full description in Ref. [25| , we can count 
the gapless quark modes by examining the spectra of 



In the vicinity of the critical temperature, chemical po- 
tential differences between colors S/i a b can be ignored as 
stated in the previous sections. 

As for the three 4x4 determinants, we can count 
the gapless modes as follows. First of all, (Rs,Bu) and 



(Bu,Bd,Rs,Gs) quarks are gapless in the dSC and 2SC 
phases, respectively, because they do not participate in 
the pairing. Even for quarks participating in the pair- 
ing, they can produce gapless excitations near the critical 
point. As an example, consider a Fermi system composed 
of two species having different Fermi momenta, Pp and 

(2) 

Pp , and a superconducting gap A of pairing with each 
other. Then the quasiparticle spectra obtained after the 
diagonalization of the matrix fermion propagator can be 
gapless when the gap is smaller than half of the difference 
of the Fermi momenta, 

\p$-pf\>2&. (93) 

In Ref. |2(j the detailed analysis of the gapless super- 
conductivity at finite temperature has been done; the 
chemical potentials and gaps have been calculated self- 
consistently at finite temperature. 

In our case, the difference of the chemical potentials de- 
pends on temperature only in higher orders so that the 
left hand side of Eq. (|93[) can be regarded as constant in 
temperature. Then, we can solve Eq. (|93|) analytically us- 
ing Eqs. J72J), (EU, and (JSOJ). When the conditions shown 
below for the mCFL, dSC and 2SC phases are satisfied, 
we obtain two more gapless modes in each phase. The 
conditions are 

• for the mCFL phase, A2 < 

• for the dSC phase, Ai < /i e . 

• for the 2SC phase, A3 < ^^e- 

For example, let us consider the temperature satisfying 
T c m < T < T c n (the dSC region). We basically have two 
unpaired quark modes (Bu and Rs). Additionally, if the 
temperature satisfies Ai < /x e , then we obtain two more 
gapless modes. 

We move on to the 6x6 determinant. In this sector, 
there are basically no unpaired quark modes in the dSC 
and mCFL phases and one unpaired quark mode (Bs) in 
the 2SC phase. In order to find additional gapless modes, 
we examine the solutions which satisfy the equation 

= ^Gd,Bs(^ = 0,q) 1/2 

= (|q|+^) 2 (|q|+^) 2 (|q|+^) 2 

[{x(x - fj, e )(x + fl e ) + xAf +(X- He)^l 

+ (.x + / i e )A2} 2 +4A?A^]. (94) 

Here x = |q| — p,, where /2 = /i — 2/i e /3 is the chemical 
potential averaged over Ru, Gd, and Bs quarks. Then 
we find 

• for the mCFL phase, no more gapless modes. 

• for the dSC phase, one more gapless mode cor- 
responding to the mixture of Ru, Gd, and Bs 
quarks. 
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FIG. 6: The number of gapless quark modes as a function of 
temperature, together with the gaps and the electron chemical 
potential. 



• for the 2SC phase, if the gap satisfies A3 < 

we obtain two more gapless modes corresponding 
to the mixtures of Ru and Gd quarks. 

In Fig. El we show the number of gapless quark modes 
together with the gaps as a function of temperature. We 
can see that the gapless modes always exist just below 
the critical temperatures T*, Tj 1 , and T™. The number 
of the gapless modes decreases as the temperature lowers 
at high density. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

In this paper, we have investigated thermal phase tran- 
sitions of the color superconducting quark matter in the 
presence of finite strange quark mass and under charge 
neutrality. We constructed a general form of the GL 
potential in such a way that it fulfills symmetry con- 
straints and then obtained the weak coupling expression 
for the GL potential including the corrections from un- 
equal quark masses (m„ i—u,d,s) and unequal quark 
chemical potentials (p,i, i=u,d,s). The corrections to 
the quadratic term of the GL potential from rrii and fii 
indicate that it is the average density of states of paired 
quarks which determines the critical temperatures of dif- 
ferent pairs, as can be seen in Eqs. (|57|l and l|58|) . 
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We found that the effects of the strange quark mass 
(m s ) through the mass term in the quark propagator 
and through the electron charge chemical potential play 
a vital role near the super-normal phase boundary, while 
the color neutrality plays a minor role in contrast to the 
case at zero temperature. 

As shown in Figs. |3] and an interplay between these 
two effects of m s splits the single phase transition (CFL 
— > normal) for the u-d-s symmetric system to three suc- 
cessive second-order phase transitions (mCFL — > dSC — > 
2SC — * normal). The window of the "dSC" phase opened 
up between the 2SC and mCFL phases is a new phase 
where only the superconducting gaps with d quarks are 
non- vanishing [Til ]. 

The obtained phase structure, which is shown to be 
valid in the weak coupling (high density) region, may be 
applicable even in the strong coupling region as long as 
the successive transitions driven by the quadratic term of 
the GL potential take place in a small interval of temper- 
ature near the super-normal phase boundary. In connec- 
tion with this, it was shown recently that the dSC phase 
discussed in the present paper may be replaced by the 
uSC phase in the low density region through a "doubly 
critical" poin t on the basis of the NJL model |2^| (see 
also Ref. [221). Moreover, if we take into account the chi- 
ral condensate at low density, some interplay between the 
broken chir al sy mmetry and the color superconductivity 
is expected 28] . How to incorporate these aspects in our 
picture based on the GL potential is an interesting open 
problem. 

In the present paper, the properties of quark and gluon 
modes such as the excitation energies were also investi- 
gated. The excitation spectra of quarks indicate that 
in addition to unpaired quark modes, whose excitation 
energy is naturally zero, more than one gapless mode 



appear just below each critical temperature as summa- 
rized in Table I. Effects of the gapless modes on physical 
phenomena such as neutron star cooling is an interesting 
problem to be examined [2IJ . 

Throughout this paper, we have studied the phase 
transitions in the mean-field level. In weak coupling, 
thermal fluctuations of the gluon fields could change the 
second-order transition to the first-order one as described 
in Refs. 0> HH- As far as the order of phase transi- 
tion is concerned, the gluon fluctuations have the fol- 
lowing effects (the detailed account will be given in our 
future publication): The second order transition, mCFL 
— > dSC, remains second order. This is because all eight 
gluons are Meissner screened across T = T] and thus 
cannot induce a cubic term in the order parameter in 
the GL potential. On the other hand, the transitions, 
dSC — > 2SC and 2SC — > normal, become weak first order 
since some of the gluons become massless at T = T c I,n 
(Table I). In order to obtain a final phase diagram we 
should study the competition between the shift of the 
critical temperature discussed in the present paper and 
that coming from the fluctuation effects. 
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